THE "TWIN JET" PLANETARY NEBULA M2-9 

Mario Livio"*^ and Noam Soker^ 



ABSTRACT 

We present a model for the structure, temporal behavior, and evolutionary status of the bipolar 
nebula M2— 9. According to this model the system consists of an AGB or post-AGB star and a 
hot white dwarf companion, with an orbital period of about 120 years. The white dwarf has 
undergone a symbiotic nova eruption about 1200 years ago, followed by a supersoft x-ray source 
phase. The positional shift of the bright knots in the inner nebular lobes is explained in terms 
of a revolving ionizing source. We show that the interaction between the slow, AGB star's wind, 
and a collimated fast wind from the white dwarf clears a path for the ionizing radiation in one 
direction, while the radiation is attenuated in others. This results in the mirror-symmetric (as 
opposed to the more common point-symmetric) shift in the knots. We show that M2— 9 provides an 
important evolutionary link among planetary nebulae with binary central stars, symbiotic systems, 
and supersoft x-ray sources. 
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1. INTRODUCTION 

M2— 9 is a spectacular, well studied bipolar nebula, usually classified as a planetary nebula 
(PN). It consists of two oppositely directed, narrow, cylindrical-like bright lobes, each extending to 
an angular projected distance of 20" from the center of the nebula (Fig. la). Farther out, there 
are two "ansae," one on each side along the symmetry axis (Schwarz et al. 1997), at 60" from 
the center, moving outward at 190 km s^^ (Schwarz et al. 1997). The inner lobes on the other 
hand, expand at a much slower speed, of ~ 40 km s~^ (Schwarz et al. 1997; Solf 2000), or maybe 
even 15 km (Trammell, Goodrich &; Dinerstein 1995; hereafter TGD). The ansae have a point 
symmetrical structure (namely, are at symmetrical points when reflected about the central star). 

The nebula and its central region show many interesting properties, which are summarized in, 
e.g. Balick (1989), TGD, Hora & Latter (1994), Phillips & Cuesta (1999), Solf (2000), and Doyle 
et al. (2000; hereafter DECS). Bahck (1989), for example, finds an Ha hnewidth of 11, 000 km s'^, 
implying the presence of a wind having a velocity of ~ 5,500 km s~^ from the nucleus. Balick 
further notes that "the spectrum of M2— 9's nucleus is more similar to the slow nova RR Tel, some 
symbiotic stars, and Scyfcrt (type 1.9) galactic nuclei..." (sec also Solf 2000). 

In view of these, and other intriguing properties (see below), it it important to identify the 
nature of this nebula and to place it in a broader evolutionary context. In the present paper we 
attempt to do precisely that. We also address specifically a perplexing positional (time dependent) 
shift of the bright knots in the inner lobes, ~ 8" from the nucleus (Allen &; Swings 1972), which 
was recently studied in detail by DECS. The shifts in the northern and southern lobes are in 
the same direction, i.e. there is a mirror-symmetry rather than a point-symmetry of the bright 
knots. The picture that emerges from many studies (e.g. Allen & Swings 1972; van den Bcrgh 
1974; Goodrich 1991; Kohoutek &: Surdej 1980; and for more information see §1 of TGD) is that of 
two UV illuminating beams in a mirror symmetry, which move around the central UV source at a 
period of ~ 100 yrs. TGD conducted a detailed study of the ionization state of the shifting knots. 
One of their conclusions was that the off-knots regions are ionized by attenuated UV radiation. 
They assumed that the lobes contain rings of higher density, and concluded that: "The UV source 
illuminates a section of the lobe, causing the rings of material to glow, producing the knots and 
explaining and naturally leading to their "fixed" north-south positions. Obscuring material near 
the star keeps one side of the lobe from being illuminated by the full force of UV source." TGD 
also estimated that the column density of the obscuring material is ~ iq^^-^sa gj^-s 

Goodrich (1991; see also Scarrott, Scarrott & Wolstencroft 1993) suggested that the radiation 
that ionizes the knots comes from a companion, and not from the central star itself. DECS, on 

the other hand, attribute the ionization to a fast jet hitting the bright spots. As we shall explain 
in the Discussion (§4), we identify some difficulties in the scenario proposed by DECS. Instead, 
wc propose in the present paper a different model, which characterizes M2— 9 as a stage in the 
evolution of interacting, detached binary systems. 
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The important properties of M2— 9, as they can be deduced from the different studies are as 
follows. 

1. The bolomctric luminosity of the central star is not known. The star is heavily obscured (e.g. 
no point source is observed in the snapshot HST survey; H. Bond 2000, private communica- 
tion). A fit to the IR spectrum gives a total flux of 1.1 x 10^ L0(D/kpc)2 (S. Kwok 2000, 
private communication), but the distance is rather uncertain (Schwarz et al. 1997 suggest 
~ 650 pc). 

2. Lim &; Kwok (2000) argue that the mass loss rate from the central source is ^ 10~^ M0 yr"-*^, 
and maybe even below 10~^ Mq yr~^. 

3. The two ansae at very large distances from the lobes expand at ~ 200 km s~^ (Schwarz et al. 
1997), a high velocity compared with that of the lobes which expand at 15-40 km s~^ (Solf 
2000). 

4. The ionizing source is very probably a WD companion (e.g. Goodrich 1991; Schwarz et al. 
1997; Lim & Kwok 2000). 

5. Material moving at ^ 5,000 km s^^ is found near the center of the nebula (Balick 1989). 
This wind was suggested to result from an accretion disk around the WD companion (e.g. 
Torres-Peimbert &; Arrieta 1998). This suggestion is consistent with the general finding that 
outflow speeds are always of the order of the Keplerian speed around the central object (Livio 
1999, 2000). 

6. The ionizing radiation escapes mainly in two directions, north and south, with an inclination 
of ~ 30° to the symmetry axis. The inclination to the symmetry axis is in the same direction 
for both "beams" (west for both, or east for both; Fig. la). The beams rotate around the 
central star with a period of 120 yrs (DECS). 

7. The ionizing radiation does reach other regions around the lobes, but it is attenuated in those 
directions (TGD). 

8. The Galactic latitude of M2— 9 is the largest in the list of 43 bipolar PNe of Corradi & Schwarz 
(1995). For an assumed distance of 650 pc it implies a Galactic height of 200 pc. This may 
imply that the progenitor of M2— 9 was not very massive compared to other bipolar PNe 
progenitors (perhaps a late A or F rather than a B star). 

We argue that all of these properties can be adequately explained in the context of a model 
in which a central binary system that is a supersoft X-ray source (SXS), is in the transition phase 
from a symbiotic nebula (possibly generated by a symbiotic nova; Solf 2000) to a PN. In particular, 
we demonstrate that a system consisting of an asymptotic giant branch (AGB) or a post-AGB star 
(with its associated slow wind) and a hot WD (which has a fast wind), in an orbital period of 
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~ 120 yrs, can ionize the "walls" of the lobes in one direction, while the ionization is substantially 
attenuated in other directions. The revolving ionizing source is explained by a combination of the 
attenuating effects of the post-AGB star's wind and by the effects of a coUimatcd fast wind (CFW) 
blown by the WD companion, which is bent by the slower wind from the post-AGB star. This type 
of flow, where the CFW is being strongly pushed upon and bent by the AGE (or post-AGB) wind 
is termed a 'weak-CFW' flow by Soker & Rappaport (2000, hereafter SROO). The bent "jets" can 
both compress the ring in each lobe, and clear a way for the ionizing radiation in the two specific 
directions around the symmetry axis. The circumbinary material has special properties (in terms 
of its optical depth, sec below) for a relatively short period of time (of the order of a few hundred 
years), and this enables the formation of the circulating ionizing beams. 

Our model provides an important link between planetary nebulae and symbiotic nebulae. 
The fact that these two types of nebulae may be related in their shaping mechanisms (and that 
the nebulae are sometimes misclassified) has been known for some time (e.g. Morris 1987, 1990; 
Corradi 1995; Corradi h Schwarz 1995; Schwarz h Corradi 1992; Corradi et al. 2000). Furthermore, 
our proposed model for M2— 9 emphasizes the important role played by binaries in the shaping of 
nebulae in general, and in the formation of jets in particular. The formation of collimated jets in 
some of these systems, as evidenced, for example, by the recent observations of He2— 90 (Sahai & 
Nyman 2000), provides important clues for the processes by which astrophysical jets are accelerated 
and collimated in general (Livio 2000). 

Our dynamical model is presented in §2, the properties of the ionizing beams are discussed in 
§3, and a summary and discussion follow. 



2. AN INTERACTING WINDS MODEL 



We consider the following configuration for M2— 9, which qualitatively describes D-type sym- 
biotic systems. D-type symbiotic systems are wide binaries (orbital periods ^ 100 yrs) in which a 
white dwarf accretes from the wind of an AGB star or a Mira variable. We propose that in M2— 9 
a WD companion accretes mass from the wind of an AGB or post-AGB star. An accretion disk is 
formed, and this disk powers a mildly collimated fast wind (CFW), or a jet (see e.g. Livio 1999, 
2000 for a review of jet formation). The general geometry is presented schematically in Figure lb. 
For a typical post-AGB star of 0.6 Mq and a WD companion of 0.8 Mq, and an orbital period of 
120 years (DBCS), the orbital separation is 27 AU. We assume a circular orbit, since an eccentric 
orbit would have caused the nebula to possess a strong departure from axisymmetry (Soker, Rap- 
paport & Kiger 2001). For a slow wind velocity of 15 km the Bondi-Hoyle accretion rate is 
~ 0.01 of the mass loss rate by the post-AGB star. The latter is ~ 10"^ yr~^ in the case of 
M2— 9 (Lim &: Kwok 2000). assuming a spherical, slow wind. The condition for the formation of 
an accretion disk around the WD can be expressed in the form ^ RwD, where Rd is the outer 
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radius of the disk that would form given the wind's specific angular momentum. The disk radius 
is given approximately by (e.g. Livio Sz Warner 1984) 

- W'W (^] ( " y ( y cm (1) 

"^'-'^ U Mo Alois cmj j bokms-V ' 

where Aftot is the total system mass, a is the separation, T4ei is the relative velocity between the WD 
and the wind, and 77 is a dimensionless parameter determined by the existing density or velocity 
gradients in the wind. Typically 77 ^ 0.25 (e.g. Ruffert 1999), which for the assumed parameters 
gives Rd/RwD ~ 25. Thus, for the assumed parameters a disk can form, but the collimation process 
may not be extremely efficient, since efficient collimation requires Ra/RwD ^ 1 (e.g. Spruit 1996; 
Livio 1999). If we assume collimation to take place by a poloidal magnetic field [of the form 
Bz ^ ('"^/^WD "f" 1)"^''^] threading the disk, then the minimum jet opening angle is given by 
6 ^ {RwD/Rdy^'^ ^ 10°. This conclusion is quite consistent with the geometry of the inner lobes 
in M2— 9, which gives an opening angle (or collimation angle, defined from the symmetry axis) of 
the CFW^ of 6c — 35°. We assume that roughly 1% of the accreted mass is blown into the CFW 
(this is consistent with low-efficiency disk outflows; Livio 1999). For the CFW velocity we take 
Vf = 5, 000 km s~i (of the order of the Keplerian velocity near the WD surface and consistent with 
observations; e.g. Livio 2000, Balick 1989). We find from these values, and from the assumption 
that the mass loss rate from the AGB, or post AGB, star in the equatorial plane is somewhat 
higher than that along the polar directions (SROO) that the ratio of the slow wind momentum flux 
(pv^) to that of the CFW is ~ 5, hence the CFW will be bent and the slow wind will flow along 
and near the symmetry axis. For the assumed parameters, the slow wind impacts the CFW at a 
large angle, and the bending is less efficient than in the case a ^ h (where h is the height of the 
impact point where substantial bending occurs; Fig. lb). Bending to the asymptotic angle 9^ will 
therefore occur at a height h = a. few x/iq ~ 100 AU where ho ~ 21 AU is given by eq. 12 of SROO 
(see Fig. lb). 

The CFW passes through a strong shock and forms a hot bubble which compresses and ac- 
celerates the slow wind with which it collides, thus forming a dense shell (e.g. SROO, Fig. lb). It 
is important to note though that by "shell" we do not refer to a closed elliptical or spherical shell, 
but rather to a "corkscrew" structure, or a "snail" type structure. Because of the relatively low 
mass loss rate from the post-AGB star, even the dense shells are rather tenuous, and hence, too 
faint to be clearly observed. If we assume that all the kinetic energy of the CFW is transferred to 
the slow wind material along the lobes, which for our assumed collimation angle is of the order of 
e ~ 0.2 of the slow wind mass, we find the shell velocity to be 



MfV^f 
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' - 190 i ^^"l ( , ) km S-. (2) 
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However, the acceleration of a dense shell by a more dilute hot bubble is prone to Rayleigh- 
Taylor instabilities, and the shell most likely will break up into many clumps. This has three impor- 
tant effects. First, the acceleration is much less efficient since the hot, shocked CFW leaks between 
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these clumps. The terminal shell (clump) velocity may thus be as low as ^^ciump ~ 50 km s"^. 
Second, the clumpy shells (a shell each orbital period) allow more of the ionizing radiation to reach 
the lobes. Finally, the leakage of the CFW keeps the directionaUty of the expanding bubble to be 
mainly radial. The clumps reach the bright rings on the lobes (where the N2 and S2 knots are 
located; see TGD) on a timescale of 

tc = 500 ( ( ^:^^?^) ~' yrs, (3) 



0.025 pcy V50 km s 

where Df is the distance to the bright rings in the lobes. The dynamical age of M2— 9 is 1, 200 yrs 
(Schwarz et al. 1997), hence the first shell reached the ring ~ 700 years ago. For an orbital period 
of 120 years this means that ~ 6 shells have reached the rings. We propose that the high density 
material in the rings was formed by these clumps (see also TGD). 

The following thing should be noted in regard to the mass loss from the WD. As the WD 
accretes from the wind, it could undergo periodic symbiotic nova outbursts (e.g. Livio, Prialnik 
&; Regev 1989). Phases of a strong WD wind could follow such outbursts (e.g. Kovetz 1998 and 
references therein). In fact, all recorded symbiotic nova eruptions were followed by an intensive 
Wolf-Rayet type wind from the hot companion (Mikolajewska 2000). The recurrence time for such 
outbursts is given approximately by (Livio 1994) 

o ^7 1 n3 I -^acc \ 



1.54 I ^] - 2.0 (^V' + 0.65 



Mq J \ Mq J V Mq 



0.7 

yrs . (4) 



Thus, for a 0.8 Mq WD we can expect strong WD wind phases roughly every ~ 6700 yrs. We 
propose that such an outburst formed the "ansae" of M2— 9 ~ 1200 years ago. The long time that 
has elapsed since the eruption implies, in the case of M2— 9, a current lower mass loss rate and a 
faster velocity than in typical post eruption symbiotics. We note that the formation of jets or any 
ballistic outflows, in short bursts, leads to a linear relation between the distance from the center 
and the outflow velocity. Such a linear relation is observed in many bipolar PNe, suggesting that 
bursts are a common phenomenon in bipolar PNe. 

Another important aspect of symbiotic-nova outbursts is the fact that following eruptions, 
the WD enters a supersoft x-ray phase, in which it is characterized by an effective temperature 
of 150,000-600,000 K (e.g. DiStefano k Rappaport 1994; Yungelson et al. 1996). During this 
phase, the ionization balance of the circumbinary material is determined by the radiation from 
this embedded supersoft x-ray source (SXS). Fast collimated outflows have been observed in three 
of the SXSs: RX J0513.9-6951 (3800 km s'^, Southwell et al. 1996; Southwell, Livio, & Pringle 
1997); RX J0019.8+2156 (~ 1000 km s-\ Becker et al. 1998); and RX J0925.7-4758 (5200 km s-^ 
Motch 1998). The mass transfer rates in these systems are not known precisely, but are probably 
above 10~^ Mq yr~^ in all three sources. 
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In principle, the WD can also be a more permanent SXS, similar in nature to systems like 
CAL 83 and CAL 87. In these systems, the accretion rate onto the WD is sufficiently high for the 
hydrogen to burn steadily (e.g. van den Hcuvcl ct al. 1992). For a WD of 0.6-0.8 Mq, the accretion 
rate for steady burning is of the order of (a few) xlO~* Mq yr~^ (e.g. Nomoto 1982), definitely 
achievable by Bondi-Hoyle accretion from the post AGB star's wind. Steady burning will result in 
a soft x-ray luminosity of 

- 1^10-6 Me yr-i J V Me J ^27 Auj [l5 km s'' J "'^ ' ' 
and an effective temperature of 

Tsxs=.5xlO^( ^^3,J^^^_. ) K. (6) 



3. THE IONIZING "BEAMS" 



In order for the soft x-rays to reach the outer visible parts of the M2— 9 nebula (and cause 
them to fluoresce), they must ionize their own path through the stellar wind of the AGB star and 
other intervening matter, and/or have a path cleared by the action of the CFW coming from the 
white dwarf. Given that much of the outer nebula apparently receives a substantial ionizing flux, 
we conclude that the radiation is able to ionize its way through the stellar wind in many if not 
most directions. A number of works have addressed the question of the ionization structure and 
concomitant line emission from an ionizing source embedded in a stellar wind (see, e.g. Nussbaumer 
& Vogel 1987; Seaquist, Taylor & Button 1984). Here we are mainly concerned with the fraction of 
the radiation flux, as a function of direction, that can pass through the wind and which is available 
to illuminate the outer nebula. This radiation flux can be estimated as follows. 

Even though the ionizing radiation is emitted by the WD companion, while the slow wind 
is blown by the AGB or post-AGB star, the geometry can still be safely assumed to be nearly 
spherical. This is a consequence of the fact that the CFW clears the slow wind close to the WD, 
up to a distance of the order of the bending distance (see previous section), of rmin ~ 100 AU, 
which is considerably larger than the orbital separation of a ~ 27 AU. For the ionizing radiation 
to be able to reach the lobes, it needs to ionize both the newly ejected slow wind material, and 
the recombining wind, in its way. For the parameters of the flow discussed here the newly ejected 
material is not important. The total recombination rate from a minimum distance rmin to the lobes 
at ri ^ Tmin IS given by (assuming a spherical flow) 

R = 4n r ccn^nydr = 8 x 10^^ ( j^f^ , ] ' ( , ) ' ( s'^ , (7) 

ir^i^ \^10-6 Mo yr-iy Vl5kms-V VlOO AUy 
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where a is the recombination coefficient, rXg and are tfie electron and ion number densities 
(respectively). Here we assumed a fully ionized nebula of solar composition, and for the density 
we used p = AfAGB[4vrr^^;s]^^. The number of ionizing photons emitted by the WD per unit time, 
assuming a temperature of Tgxs — 5 x 10^ K, is 

S ~ 4 X 10^6 ( \ (8) 

Vl037 erg s-V ^ ^ 

The two rates obtained above, for R and S, are within an order of magnitude from each 
other. This implies that a substantial fraction, but not all, of the radiation will be absorbed along 
directions facing the AGB star, while only a negligible fraction will be absorbed along directions 
cleared by the CFW. The CFW clears away the slow wind along its path. If the compressed shell 
moves at ~ 50 km s^^ (eq. 2 and following discussion), then during about half an orbital period 
(60 years), before fresh slow wind material enters the swept-up region, it will move to a distance 
of r ~ 600 AU. This means that the total number of recombinations along that direction is much 
smaller than the number of ionizing photons emitted in that direction. Along other directions, the 
slow wind's closest distance to the WD will be less than the distance of bending, and a substantial 
fraction of the ionizing flux will be absorbed, but some will still reach the lobes. This is consistent 
with deductions from observations (TGD). The estimates presented by cqs. (7) and (8) therefore 
support our proposed scenario, in which the east-west motion of the bright knots is caused by 
ionizing beams which pass through channels cleared by bent "jets" (or a CFW). 



4. DISCUSSION AND SUMMARY 



M2— 9 is a unique nebula in that it shows a variability on a time scale which is consistent 
with the orbital period expected from the binary model for the formation of bipolar nebulae. In 
this model most (although not all) bipolar nebulae axe formed by binary systems with orbital 
periods in the range of ~ 1-10^ yrs (e.g. SROO). The time varying features exhibit a mirror 
symmetry, rather than a point-symmetry (TGD; DECS). As we have shown, a mirror symmetry 
is expected from the interaction of two winds, a slow wind emitted by the AGB or post AGB 
star, and a jet (or collimated fast wind [CFW]) emanating from the accretion disk around the WD 
companion. The mirror symmetry is produced as the CFW is pushed and bent by the slow wind. 
Our model identifies the evolutionary state of M2— 9 as one which relates three different classes of 
binary systems: supersoft X-ray sources (SXSs), symbiotic nebulae, and PNe (or proto-PNe). The 
ionizing source is the white dwarf companion, and not the mass-losing star (as already concluded 
by e.g. Goodrich 1991; Schwarz ct al. 1997; Solf 2000; Lim & Kwok 2000). At the same time, the 
fast outward moving "ansae" (Schwarz et al. 1997) and the large and bright nebula strongly suggest 
that in the recent past the system has both undergone a more intense mass loss episode, and an 
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outburst associated with the accreting WD companion. Hence, our evolutionary scenario for the 
nebula can be characterized as a symbiotic nebula evolving to a PN. 

While M2— 9 is unique in its time variation, it is not in its morphology. The bipolar proto 

PN Henize 401, for example, has a similar structure (Sahai, Bujarrabal, & Zijlstra 1999), but with 
no evidence for circulating ionizing beams. Wc note that Sahai et al. (1999) argue that the bipolar 
structure in He 401 was formed by collimated bipolar jets, and not by the more commonly assumed 
interacting winds model, in which a spherical fast wind catches up with a previously ejected bipolar 
slow wind (e.g. Balick 1987). The model of Sahai et al. has in fact some features that are similar 
to those of our model for M2— 9. 

Our model offers an interpretation to the bright knots that is different from the one proposed by 
DECS. We showed that the bright knots are exposed to a higher flux of ionizing radiation (from the 
hot WD) than the rest of the lobes, by virtue of the path cleared by the CFW. DECS, on the other 
hand, attribute the brightness/ionization of the knots to shock heating, caused by jets impacting 
on lobe walls. We see a few problems with this suggestion: (i) it offers no clear explanation for the 
mirror symmetry (the jets will collide with slow wind material ejected half a period earlier and will 
be slowed down substantially), (ii) The fact that the off-knots regions are ionized as well (TGD) 
argues for an ionizing source that is not confirmed to the small cross section of a jet, which is present 
at all times, (iii) We find the 'corkscrew' model of DECS (see e.g. their Fig. 2) to provide for a 
rather poor description of the nebula (e.g. the lines that are supposed to represent the emission in 
1989 and 1999 do not follow very well the actual emission). We rather interpret the 1989 and 1999 
images as follows. The bright emission regions (from the center of the nebula to the bright blobs 
along the symmetry axis) are moving in the same sense, with the 1989 emission being toward the 
outskirts (east side) relative to the 1999 emission. The change in the east-west location occurs only 
near the bright blobs (N3 and S3), and even there it represents a complicated flow structure, and 
radiative transfer from these regions outward. Nevertheless, our model does share some common 
features with that of DECS (e.g. a fast outflow from the companion). 

One point in our model deserves further discussion. A cleared path in a certain direction for 
the ionizing beam (and therefore enhanced ionizing radiation), could result, in principle, in a more 
extended, rather than brighter, nebula (which is not observed). We should therefore check that for 
the parameters of M2— 9 there is no significant penetration of the ionization front in the direction 
perpendicular to the lobe walls. For an electron density in the lobes of rie ~ 10^ cm~^ (DECS), an 
ionizing source of S* ~ 4 x 10^^ (eq. 8), and a path clearing for about 10% of the period, one can 
easily show that the penetration of the ionization perpendicular to the wall would be by less than 
~ O.ID (where D is the radial distance to the S2, N2 knots). Thus, while we do not expect any 
significant perpendicular penetration, we do expect the bright regions to extend along the lobes 
(due to the lower density), as observed (Fig. la). 
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To conclude, the spectacular nebula M2— 9 provides an opportunity to identify an evolutionary 
link among three types of (empirically defined) binaries: PNe with binary central stars, symbiotic 
binaries, and supersoft x-ray sources. 
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FIGURE CAPTIONS 
Fig. la: The general structure of M2— 9. 

Fig. lb: A schematic drawing of the flow structure (not to scale). The figure shows the plane 
perpendicular to the equatorial plane, and momentarily containing the two stars. The 3D structure 
is that of a corkscrew. Drawn on the figure are: the slow wind from the post-AGB star (stream 
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lines originating from the post- AGB star) ; stream lines of the pre-shocked CFW; the shock front of 

the CFW; The hot bubble region (again, a corkscrew structure in 3D); and the dense shell formed 
from the slow wind material accelerated by the hot bubble, with its Rayleigh- Taylor instabilities. 
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